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A global analysis of the experimental data on azimuthal asymmetries in semi-inclusive deep inelas- 
tic scattering (SIDIS), from the HERMES and COMPASS Collaborations, and in e + e" -> hxhiX 
processes, from the Belle Collaboration, is performed. It results in the extraction of the Collins 
fragmentation function and, for the first time, of the transversity distribution function for u and 
d quarks. These turn out to have opposite signs and to be sizably smaller than their positivity 
bounds. Predictions for the azimuthal asymmetry A!^f;^ s+ ^ h \ as will soon be measured at JLab 
and COMPASS operating on a transversely polarized proton target, are then presented. 



PACS numbers: 13.88.+e, 13.60.-r, 13. 66. Be, 13.85.Ni 
I. INTRODUCTION 

The transversity distribution function, usually denoted 
as hi q (x,Q 2 ) or Atq(x,Q 2 ), together with the unpolar- 
ized distribution functions q(x,Q 2 ) and the helicity dis- 
tributions Aq(x, Q 2 ); contains basic and necessary infor- 
mation for a full understanding of the quark structure, 
in the collinear, k± integrated configuration, of a polar- 
ized nucleon. The distribution of transversely polarized 
quarks in a transversely polarized nucleon, ATq(x,Q 2 ), 
is so far unmeasured. The reason is that, being related 
to the expectation value of a chiral-odd quark operator, 
it appears in physical processes which require a quark 
helicity flip: this cannot be achieved in the usual inclu- 
sive DIS, due to the helicity conservation of perturbative 
QED and QCD processes. 

The problem of measuring the transversity distribu- 
tion has been largely discussed in the literature The 
most promising approach is considered the double trans- 
verse spin asymmetry Att in Drell-Yan processes in pp 
interactions at a squared cm. energy of the order of 
200 GeV 2 , which has been proposed by the PAX Col- 
laboration @, H, H, @. However, this requires the avail- 
ability of polarized antiprotons, which is an interesting, 
but formidable task in itself. Meanwhile, the most ac- 
cessible channel, which involves the convolution of the 
transversity distribution with the Collins fragmentation 
function [||, is the azimuthal asymmetry A s ^^ s+ ^ h ^ 1 in 
SIDIS processes, namely £p^ —> tirX. This is the strat- 
egy being pursued by HERMES, COMPASS and JLab 
Collaborations. 

A crucial improvement, towards the success of this 
strategy, has been recently achieved thanks to the inde- 
pendent measurement of the Collins function (or rather, 
of the convolution of two Collins functions), in e + e~ —> 
h\hiX unpolarizcd processes by Belle Collaboration at 
KEK [7]. By combining the SIDIS experimental data 



from HERMES [|©] and COMPASS [3, with the Belle 
data, we have, for the first time, a large enough set of data 
points as to attempt a global fit which involves, as un- 
known functions, both the transversity distributions and 
the Collins fragmentation functions of u and d quarks. 

In Section |TT] we briefly remind the basic formalism 
involved in the description of the SIDIS asymmetry 
^sm{ct>s+4>h) ^ anc j j n s cc ti on [nil we develop, in somewhat 
greater detail, a similar formalism for the azimuthal cor- 
relations, involving two Collins functions, measured by 
Belle in e + e~ — > hxh^X processes. In Section ITVl we 
perform a slobal fit of HERMES COMPASS [TcJ 

and Belle [7J data, in order to extract simultaneously the 
Collins fragmentation function A N D 7r / q y (z,p±) and the 
transversity distribution function Axq{x) for q = u,d. 
We then use, in Section [Vj the transversity distributions 
and the Collins functions so determined, to give predic- 
tions for forthcoming experiments at JLab and CERN- 
COMPASS. Comments and conclusions are gathered in 
Section EI 



II. TRANSVERSITY AND COLLINS 
FUNCTIONS FROM SIDIS PROCESSES 

The exact kinematics for SIDIS £p — ► I h X processes 
in the 7* — p cm. frame, including all intrinsic motions, 
was extensively discussed in Ref. and is schemati- 
cally represented in Fig. [TJ We take the virtual photon 
and the proton colliding along the z-axis with momenta 
q and P respectively, and the leptonic plane to coincide 
with the xz plane. We work in the kinematic regime in 
which Pt — Aqcd — k±, where k± is the magnitude of 
the intrinsic transverse momentum k± of the initial quark 
with respect to the parent proton and Pt = \Pt\ is the 
magnitude of the final hadron transverse momentum. We 
neglect second order corrections in the k±/Q expansion: 
in this approximation, the transverse momentum of 
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FIG. 1: Three dimensional kinematics of the SIDIS process, 
according to Trento conventions [l4| ■ The photon and the pro- 
ton collide along the i-axis, while the leptonic plane defines 
the xz plane. The fragmenting quark and the final hadron 
h are emitted at azimuthal angles <p and 4>h, and the proton 
transverse spin direction is identified by <f>s- 



the observed hadron h with respect to the direction of 
the fragmenting quark is related to k± and Pt by the 
simple expression p ± = PT — zk±; in addition, the light- 
cone momentum fractions x and z coincide with the usual 
measurable SIDIS variables, z = Zh = {P ■ Ph) / [P ■ q) and 
x = Xa = Q 2 /(2P ■ q). In this region factorization holds 
(l2l . Il3j , leading order iq^> Iq elementary processes are 
dominating and the soft Pt of the detected hadron is 
mainly originating from intrinsic motions. 



The transverse single spin asymmetry (SSA) for this 
process is defined as 



d 6 c 



J^t'hX 



r i P ^i'hx d(jT _ da i 



*fhX + d 6 c 



where d 6 a e 



l -*£hX — 



(1) 

dcr'^ is a short hand notation for 
(d 6 a l P r ' l ^ hx )/(dx B dydz h d 2 P T d<f) S ). It will often hap- 
pen, in comparing with data or giving measurable pre- 
dictions, that the numerator and denominator of Eq. ((T|) 
will be integrated over some of the variables, according 
to the kinematical coverage of the experiments, f and 
I refer, respectively, to polarization vectors S and —S, 
see Fig. Q] A full study of Eq. (Q]) , with all contributions 
at all orders in k±/Q, will be presented in a forthcoming 
paper [l5j . 

We consider here, at 0(k±/Q), the sin(0s + 4>h) 
weighted asymmetry, 



^jsin(0 s + 0h) 



[dcr T - da 1 ] sin(0 s + <j> h ) 



J d<f>s d(f)h [rfcr T + do 4 ] 



measured by the HERMES @, ©] and COMPASS [nj 
Collaborations. This asymmetry singles out the spin de- 
pendent part of the fragmentation function of a trans- 
versely polarized quark with spin polarization s and 
three-momentum p q : 

Dh/q,s(z,P ± ) = D h/q (z, P± ) 

+ -A N D h/q t(z, Px )s-(p q xp 1 _), (3) 

resulting in 



J 



2 elf, #s # h d 2 k ± A T q(x, k±) ^ N D h/q r (z,pi_) sin(0 s +(p + $) sin(0 5 + cf> h ) 



Asin(<t>s+(j> h ) _ 1 

UT 



f da 
^2 e l J d<j) S d(j) h d 2 k ± _ f q /p(x,k±) — D h / q {z,pj_) 

I 



(4) 



In the above equation ATq(x,k±) is the unintegrated da/dy is the planar unpolarizcd elementary cross section 
transvcrsity distribution, 



A T q{x) = h lq {x) = / d z k ± A T q(x, k x ) , (5) 



while A Dji/qi (z,p±) is the Collins function, often de- 
noted as LL 



A N D h/q r(z, Pl _) = ^ H^(z, Px ) 
zrrih 



(6) 



da 2ira 



dy sxy 1 



[l + (l-y) 2 ] 



(7) 



and 

d(Aa) da 1 "' 



da^ 



dy 



dy 



dy 



47ra 2 ,„ 

— 2 i-y • (8) 

sxy 



The sin (05 + ip + <fi q ) azimuthal dependence in Eq. (H} 
arises from the combination of the phase factors in 
the transvcrsity distribution function, in the non-planar 
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£ q — > I q elementary scattering amplitudes, and in the 
Collins fragmentation function; <j)g and ip identify the di- 
rections of the proton spin S and of the quark intrinsic 
transverse momentum k±, see Fig.[T] (j> q is the azimuthal 
angle of the final hadron h, as defined in the fragmenting 
quark helicity frame. Neglecting 0(k 2 _/Q 2 ) terms, one 
finds 



cos< 



sin 



COS(0fc, —<p)—z — , 

p± p± 



Pi 



sin((f>h - (p) 



(9) 



A full study of Eq. ([2]) , taking into account intrinsic mo- 
tions with all contributions at all orders, following the 
general approach of Ref . [l6| , will be presented in a forth- 
coming paper (l5| . Here, in agreement with all papers on 
the Collins effect in SIDIS so far appeared in the litera- 
ture, we work at 0(k±/Q) and use Eqs. ((4|) and ([9]). 

fq/ p (x,k±) is the unpolarized transverse momentum 
dependent (TMD) distribution function of a quark q in- 
side the parent proton p, while Dh/ q {z,p±_) is the unpo- 
larized TMD fragmentation function of quark q into the 
final hadron h. We assume the k± and p± dependences 
of these functions to be factorized in a Gaussian form, 
suitable to describe non-perturbative effects at small Pt 
values and simple enough to allow analytical integration 
over the intrinsic transverse momenta: 



f q / p (x,kj_) = f q /p{x) 



D h/q( Z ,P±-) = P>h/q(z) 



-pI/(pI) 



(10) 
(11) 



where f q / p (x) and Df l / q {z) arc the usual integrated par- 
ton distribution and fragmentation functions, available in 
the literature; in particular we refer to Refs. [13, HU and 
fl9| . The QCD induced Q 2 dependence of these functions 
is also taken into account, although we do not indicate 
it explicitly. Finally, the average values of k\ and p\ 
arc taken from Ref. pd| . where they were obtained by 
fitting the azimuthal dependence of SIDIS unpolarized 
cross section: 

{k\) = 0.25GcV 2 , (pi) = 0.20 GcV 2 . (12) 

Notice that such values are assumed to be constant and 
flavor independent. 



The transversity distributions and the Collins func- 
tions are unknown. We choose the following simple pa- 
rameterization 



-k 2 ± /(kl) q 



A T q(x, k X ) = - Mq (x) [f q/p (x) + Aq(x)] , 



(13) 



A N D h/q i(z, P± ) = iMq (z) D h/q {z) h(p±) 



(14) 



with 



Nq (x) = N? x a (l - xf , (15) 



NZ(z)=N§ z<(l- z) s U ^°> 



1 



i8 s 



(16) 
(17) 



and |7V«j r |, |iV g c | < 1. In general (k\) T ^ (k 2 ± ), but from 
our fits we learn that present experimental data are in- 
sensitive to such a difference, therefore we simply assume 
(kj_) T = (kj_). Also, in this first simultaneous extraction 
of the transversity and Collins functions, we let the co- 
efficients N q and N q to be flavor dependent (q — u,d), 
while all the exponents a, (3, 7, S and the dimensional pa- 
rameter M are taken to be flavor independent. 

Notice that our parameterizations are devised in such 
a way that the transversity distribution function auto- 
matically obeys the Soffer bound [2(| 



\A T q(x)\ < - [f q/p (x) + Aq(x)} , (18) 



and the Collins function satisfies the positivity bound 

\A N D h/q i (z,px)\ < 2D h/q {z,p±) , (19) 

since J\f q (x), Nq (z) and h(p±) are normalized to be 
smaller than 1 in size for any value of x, z and p± re- 
spectively. 

By insertion of the above expressions into Eq. ([4]), we 
obtain, in agreement with Refs. pll. I22I], 



Pt 1-y /o- c (P±Y: 



UT 



M sxy 



e 2 q 



(Pt)1 



1 



(x) [f q/p (x) + Aq(x)] (z) D h/q {z) 



(P T ) sxy 2 



e-lf q / P (x) D h/q {z) 



(20) 
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where 

(21) 

Eq. expresses ,4™ ) m terms of the parame- 

ters a, (3, 7, <5, iVg , Nq and M. In Section [W] we shall fix 
them by performing a best fit of the measurements of 
HERMES, COMPASS and Belle Collaborations. Actu- 
ally, we shall consider, following the experimental data, 
^jsmOs+<M as a f unc tion of one variable at a time, by 
properly integrating the numerator and denominator of 
Eq. (|20|) : the integration over x and z gives the Pt distri- 
bution of A s "^ s+ ^ h ' > , whereas the integrations over Pt 
and z or Pt and x, yield the x and z distributions. Notice 
that, with our approximations, x = x B and z = Zh- 




FIG. 2: Three dimensional kinematics of the e e~ — > hi/12 X 
process, in the gg cm. frame. In this configuration the re- 
constructed thrust axis identifies the ^-direction, the lepton- 
quark scattering plane defines the xz plane. 



III. COLLINS FUNCTIONS FROM e+e" 
PROCESSES 

The kinematics corresponding to the e + e~ — * h\h 2 X 
process is schematically represented in Fig. [5] the two 
detected hadrons h\ and h 2 are the fragmentation prod- 
ucts of a quark and an antiquark originating from e + e _ 
collisions. We choose the reference frame so that the 
e + e~ — > 99 scattering occurs in the xz plane, with the 
back-to-back quark and antiquark moving along the z- 
axis. This choice requires, experimentally, the recon- 
struction of the jet thrust axis, but it involves a very 
simple kinematics and a direct contribution of the Collins 
functions, as we shall see. A different choice, originally 
suggested in the literature [23j], is discussed at the end 
of this Section. In the configuration of Fig. [21 the four- 
momenta of the e + ,e~ (fc + ,fc~) and of the q, 9 (91,92) 
are 



91 = ^(1,0,0,1), 9 2 = ^(1,0,0,-1), (22) 

fc~= -^(1 -sin 0,0, cos 0), fe+=^-(l,sin0,O-cos0) . 

The final hadrons h\ and hi carry lightcone momentum 
fractions z\ and zi and have intrinsic transverse momenta 
p^j and p^ 2 with respect to the direction of fragmenting 
quarks, 

P_li = pj.i(cosvi,sinvi,0), 

P±2 = p±i(cosipi,sinip 2 ,0) , (23) 

so that their four-momenta can be expressed as 

Pi = [zi — + , p±i cos ip 1 , p ± i sinip 1 , 

V 2 2ziVs 



2 2z ly /ZJ ' 



(24) 



P2 = ( Zi — ^ + n P±2 r - , P±i cos 952 , p±i sin ip 2 , 
V 2 2z 2 Js 



-z 2 



P±2 
2z 2 y r s 



(25) 



At large cm. energies and not too small values of z, 
one can neglect second order corrections in the p±/ (zyfs) 
expansion, to work with the much simpler kinematics: 



V s • V s 

-Pi = \Zl — , pj_iCOsyi, p_Li sm^i, Zl — 



(26) 



P2 = {Z2 — , p±2 COS(/3 2 , RL2 SU1(£2, ~Z2 — J • (27) 

Notice also that in this limit the lightcone momentum 
fractions z coincide with the observable energy fractions 

Zh, 



z h = 2E h /^~s =z + ?±~z 
zs 



(28) 



The cross section corresponding to this process, with 
unpolarizcd leptons, can be written as 



^.e+e ->h 1 h 2 X 



(29) 



dz\ dz 2 d 2 p ±1 d 2 p^ 2 d cos 6 

q {A} 

x< 1 ^(«i > P±i)^{f(^ > P± 2 ) > 

where Ma,a ? ;A + a_ are the helicity amplitudes corre- 
sponding to the elementary scattering e + (A+)e~(A_) — > 
9(A g )9(Ag), 9 = u,u,d,d, s,s (neglecting heavy flavors) 
and X){a} indicates a sum over all helicity indices. In 
this case there are only two non-zero, independent am- 
plitudes: 



M 4 



= M_+._+ = e 2 e g (l + cos0) , 



M-+:+- = M+_._ + = e 2 e q (l - cos (9) . (30) 

The functions D^Jy (zi , p ±1 ) and D^ 2 jj_ (z 2 ,p ±2 ) are the 
probability densities which describe the fragmentation 
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of quarks and antiquarks into the physical hadrons hi 
and h 2 respectively (see Section II. C of Ref. [l6| for de- 
tailed explanations). In particular, the diagonal elements 

-D++(z,p^) and D h / q _(z 1 p^) correspond to the trans- 
verse momentum dependent unpolarized fragmentation 
function D h / q (z,p±), 

D h + ^(z,p ± ) = D h _^(z, P± ) = D h/q (z, P± ) , (31) 
whereas the non-diagonal elements 



D h J«(z, P± ) = D'^(z, Pl _) 



h/q. 



(32) 



D^(z,p x ) = D n _'«(z, Ps _)e-^ = -Byi(z,p ± )e- i * , 

are related to the Collins fragmentation function 
A N D h/q ,(z,p x ) [IH by 



■)h/g, 



A N D h/g ,(z,p x ) = -2iD h Jl(z, Pl _) = 2iD h Jl(z, Pl _) 



(33) 

The angle tp in Eq. (|33[) is the azimuthal angle identifying 
the direction of the observed hadron h in the helicity 
frame of the fragmenting quark q. Similar relations hold 
for the antiquark fragmentation functions, where one has 
to take into account a sign difference in (p originating 
from the fact that the antiquark is chosen to move along 
the — z direction. Finally, inserting Eqs. (|3TH) - ([3"3")l into 
Eq. (|29[) and performing the sum over the quark hclicitics 
one obtains 



da 



e^e — >h\h 2 X 



dz\ dz2 d 2 p ±1 d 2 p j_ 2 dcos 6 



3na 2 
2s 

1 



^2 e q {(1 + cos 2 6) D hl/ g(zi,p_Li) D h2/ g(z 2 ,px2) 
q 

A N D hl/q7 {zx,pj_i) A N D h2 /^ {z 2 ,p± 2 ) cos(<£i + ip 2 )\. 



sin 2 



Eq. (j34|) shows that the study of the correlated produc- 
tion of two hadrons (one for each jet) in unpolarized e + e~ 
collisions offers a direct access to the Collins functions, 
both regarding their z and p± dependences. So far, only 
data on the z dependence are available. Notice that by in- 
tegrating over the intrinsic transverse momenta p_ L1 and 
Pj_ 2 one recovers the usual unpolarized cross section, 



da e 



*hxh 2 X 



dz\ dz2 d cos 

37TO 2 



(35) 



2s 



(1 + cos 2 6) 4 D hl/q ( Zl )D h2/q -{z2) , 



having used 



J d 2 p ± D h/q (z, PJ _) = D h/q (z 



(34) 



(36) 



To construct the physical observable measured by the 
Belle Collaboration, we now perform a change of angular 
variables from (951, ^2) to ((/?i, ipi+^2) arid then integrate 
over the moduli of the intrinsic transverse momenta, p±i 
and pj_2, and over the azimuthal angle (p\. This leads to 



^.e+e — >h 1 h 2 X 



dzi dz2 dcosu d((pi + ip2) 4s H I 

+ i sin 2 cos(^i + tp*) A N D hl/qT (zi) A N D h2 / qr (z 2 )} 



r 



(37) 



where we have defined 



By normalizing Eq. (|37[) to the azimuthal averaged cross 
section, 



d 2 Px A N D h/q ,(z, Pl _) ee A N D h/q r(z) 



(38) 



(da) 



]_ da ^h t h 2 x 
2tt dz\ dz2 d cos 9 
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= 17 E e \ 0- + cos2 °) D hl/q (zi) D h2/q {z 2 ) , (39) 
one has 

A(z 1 ,z 2 ,0,ip 1 + Lp 2 ) = 



1 



>h 1 h 2 X 



= 1 + T 



1 sin 



(ofcr) c?zi dz 2 dcos 9 d(<pi + (^ 2 ) 
— cos(^i + ip 2 ) 



4 1 + cos 



E q e 2 q A N D hl/q r( Zl )A N D h2/q r(z 2 ) 



(40) 



Actually, Belle data are collected over a range of 9 
values, according to the acceptance of the detector (see 



Eq. 1(63])). Thus, Eqs. ([37]) and ((39j) are integrated over 
the covered range resulting in some specific (sin 2 9) and 
(1 + cos 2 9) values. 

Finally, to eliminate false asymmetries, the Belle Col- 
laboration considers the ratio of unlikc-sign to like-sign 
pion pair production, Ajj and Al, given by 



R = 



M 

A, 



1 



+ J C0S(<^1 + ip 2 ) -JJ 



V p 
+cos 2 9) r U 



1 + i cosfai + y 2 ) (1 < +" 2 s 2 ) e > Pl 
(sin 2 L 

<P2) 



1 



1 H — cos(g9i , 

4 m r ' (1 + cos- 

1 + cos(^i + Lp 2 ) A 12 {z 1 ,z 2 ) 



(Pu-Pl) 



(41) 



with 



J 



Pn 



Pr 



E q 4 [A N D„ +/q r (zi) A"D*-,gr (z 2 ) + A N D„- /q r (zi) A N D„ +/q , (z 2 )] 

T,q e 'q [ D Tt+/q{ Z l) D v - /q (z 2 ) + D v -/ q {z{) D„+ /q (z 2 )] 

e 2 [A* (zQ A w J} T+/i -, (z 2 ) + A N D v - /q i (zQ A N D 7T - /q r (z 2 )} 



A 12 (zi,z 2 ) 



1 



T, q e l [D^+/ Q {zi) D v+/q {z 2 ) + D„- /q {zi) D n -/ q {z 2 ) 
(Pu - Pl) ■ 



4(1 + cos 2 



(42) 
(43) 
(44) 



For fitting purposes, it is convenient to re-express Pjj and Pl in terms of favoured and unfavoured fragmentation 
functions. 



7T+ ju 
7T+/d 



D 



TV+/d 



-Id 



D, 



— P > TT+/U — D-rr- /„ — D 



TT I U 



-/a 

y 7T-/d 



D 

= D 



fav 



7T±/S 



unf 3 



and similarly for the A^D, obtaining 

[5 A N D iav ( Zl ) A N D iav (z 2 ) + 7A N D unf ( Zl ) A N D unf (z 2 ) 



Pu = 
Pl = 



[5 D {sv (zi) D ia , v (z 2 ) + 7 D un f (Zl ) £) unf (z 2 )] 

[5 A N D iav ( Zl ) A N D nni (z 2 ) + 5 A N D unf ( Zl ) A N D fav (z 2 ) + 2 A N D VBt (z 1 ) A N D unf (z 2 )} 
[5 Df av (zi) D unf (z 2 ) + 5D unf (zi) -Dfav(z 2 ) + 2L» unf (zi) D unf (z 2 )] 

I 



(45) 
(46) 



(47) 
(48) 



having neglected heavy quark contributions. Pjj and Pl 
arc the same as in Ref. [13], remembering Eq. ((6]) and 
noticing that 



A N D h/ql (z)=Jd 2 p ± A N D h/q ,(z, P± ) 



J zmh 



-L(l/2)g 



(49) 
(z). 



In addition, the Belle Collaboration presents a second 
set of data, analysed in a different reference frame: fol- 
lowing Ref. [IK, one can fix the z-axis as given by the 
direction of the observed hadron h 2 and the xz plane as 
determined by the lepton and the h 2 directions. There 



will then be another relevant plane, determined by z and 
the direction of the other observed hadron hi, at an an- 
gle <pi with respect to the xz plane. This kinematical 
configuration is shown in Fig. [3J it has the advantage 
that it does not require the reconstruction of the quark 
direction. 

However, in this case the kinematics is more compli- 
cated. At first order in p±/(zy/s) one has 



P 2 = 1^1(1,0,0,-1), 



(50) 



, V* P±2 P±2 . Vs\ , K1 s 
q 2 = — , COS^ 2 , Sin^ 2 , 7T ) ( 51 ) 

2 z 2 z 2 2 J 
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/ \fs pj_2 P±2 . \AM /rr.\ Moreover, the elementary process e + e — > no docs not 

\ 2 Z2 22 2 / occur in general m the xz plane, and thus the hchc- 



/.s 



ity scattering amplitudes involve an azimuthal phase ip 2 - 
Pi = [ Pit cos fa , Pit sin fa , Z\ —— J , (53) One can still perform an exact calculation, using the gen- 

eral approach discussed in Ref. [l6j]. A detailed descrip- 

PitCOS^i - — p± 2 cos ip 2 , tion wil1 bc P rcscntcd m a forthcoming paper [l5|. We 

Z2 give here only the results valid at 0{p±/z-\/s). The ana- 

Zl \ logue of Eq. (f3"4"| now reads 

Pit sin 0! pj_ 2 smip 2 , . (54) 

22 J 



P±i 



da e+e ^hih a X 3na 2 



dz 1 d Z 2d^ 1 d^ 2 dco S 02 ~ 2s E4{( 1 + «« a ^)^/ fl (*i.J'xi)U W ,(*,PJ J ) (55) 

+ J sin 2 2 AX/^^-PiO^^^rfe.P^) cos(2p 2 + 



where fa^ 1 is the azimuthal angle of the detected hadron 
hi around the direction of the parent fragmenting quark, 
q. Technically, fa^ 1 is the azimuthal angle of p ±1 in the 
hclicity frame of q. It can bc expressed in terms of the in- 
tegration variables we are using, p ±2 an d Pit- At lowest 
order in p±/(z^/s) we have 



cos 



L 0g 1 = 



Pit 
P±i 
Pit 

P_L1 



cos 



Zi p± 2 

tP2) , 

Z2 RL1 



sm(0i - ip 2 ) 



(56) 
(57) 



Integrating Eq. (|55[) over pj_ 2 an d Pit, but not over fa, 
and normalizing to the azimuthal averaged unpolarizcd 
cross section (|39p , we obtain the analogue of Eq. (TTOl) , 




FIG. 3: Three dimensional kinematics of the e + e~ — > /u/12 X 
process. In this configuration the z direction is identified by 
the momentum of the final hadron /12, while hi is emitted at 
an azimuthal angle fa with respect to the lepton-/i2 plane, 
defined as the xz plane. 



A( Zl ,z 2 ,e 2 ,fa) = 1 + - 



sin 2 9 2 



+ zk 1 + cos 2 : 



cos(2 fa) 



E 



£A N D hl/q r( Zl )A N D h2/(j i(z 2 ) 



q Q 



E 



q elD hl/q { Zl )D h2/ij (z2) 



(58) 



in agreement with Ref. [24[ taking into account the dif- 
ferent notations, Eqs. ^ and ([35]). 

Finally, Eq. (|4"Tj) becomes in this configuration 



i?~ l + cos(2^i)^o(zi,« 2 ) 



(59) 



with 



1 z\ z 2 (sin 62) , „ 

7T Zf + Zf (1 + COS 2 6 2 ) 

where P\j and Pj, are the same as defined in Eqs. (14"7]) 
and 



IV. TRANSVERSITY AND COLLINS 
FUNCTIONS FROM A GLOBAL FIT 



We can now pursue our strategy of gathering simulta- 
neous information on the transversity distribution func- 
tion Axq(x,k±) and the Collins fragmentation function 
A N D h / q i (z,p±). To such a purpose we perform a global 
best fit analysis of experimental data involving these 
functions, namely the data from the SIDIS measurements 
by the HERMES @, @ and COMPASS Q3] Collabora- 
tions, and the data from unpolarizcd e + e~ — > ft,i/i 2 A" 
processes by the Belle Collaboration 0]. 

Axq(x, k±) and A N D h / q i (z,p±) are parameterized as 
shown in Eqs. (| 1 3 [) — (fTTj) . Considering the scarcity of 
data, in order to minimize the number of parameters, we 
assume flavor independent values of a and (3 and, simi- 
larly, we assume that 7 and 5 are the same for favored and 
unfavored Collins fragmentation functions, Eqs. (|45[) and 
f|4(3[) ; we then remain with a total number of 9 parame- 
ters. Their values, as determined through our global best 
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TABLE I: Best values of the free parameters for the u and 
d transversity distribution functions and for the favored and 
unfavored Collins fragmentation functions, Eqs. (|13|1 - (|17[) . as 
obtained by simultaneously fitting HERMES and COMPASS 
data on the A^ 't's+^h) aS y mme try and the Belle data on the 
A12 asymmetry, Eq. (|44p , proportional to cos^i+i^)- Notice 
that the errors generated by MINUIT are strongly correlated, 
and should not be taken at face value. The significant fluctu- 
ations in our results are shown by the shaded areas in Figs. [4] 
[5]and[6] as explained in the text. The values of {k 2 _) T = (k\) 
and (p 2 ± ) are fixed, according to Eq. (IT2")) . 



TABLE II: Best values of the free parameters for the u and 
d transversity distribution functions and for the favored and 
unfavored Collins fragmentation functions, Eqs. (|13|1 - (|17[) . as 
obtained by simultaneously fitting HERMES and COMPASS 
data on the < ^ s+< < ! ''») asymmetry and the Belle data on the 
Aq asymmetry, Eq. (|60|l . proportional to cos(2</>i). Notice 
that the errors generated by MINUIT are strongly correlated, 
and should not be taken at face value. The significant fluctu- 
ations in our results are shown by the shaded areas in Figs. [4] 
[5]and[6l as explained in the text. The values of (fcj_)r = {k\) 
and (p\) are fixed, according to Eq. ()12p . 



FIT I (A 12 ) 



X 2 /d.o.f. 



0.81 



Transversity N£ = 0.48 ± 0.09 N} = -0.62 ± 0.30 
distribution a =1.14 ± 0.68 /3 = 4.74 ± 5.45 
function (k 2 ± ) T = 0.25 GeV 2 



FIT II (A)) 



X 2 /d-o.f. 



0.77 



Transversity 
distribution 
function 



N 



= 0.42 ±0.09 Nj 
a = 1.20 ±0.83 P 

2 



(ki) T = 0.25 GeV 



-0.53 ± 0.28 
5.09 ± 5.87 



Collins iV£ v 
fragmentation 7 
function {p\) 



= 0.35 ±0.16 
= 1.14 ±0.38 
= 0.20 GeV 2 



Id 
S 

M 2 



-0.85 ±0.36 
0.14 ±0.36 
0.70 ±0.65 
GeV 2 



Collins JV£, 
fragmentation 7 
function (p^ 



= 0.41 ±0.10 
= 0.81 ±0.40 
= 0.20 GeV 2 



c 



5 

M 2 



-0.99 ± 1.24 
0.02 ±0.37 
0.88 ± 1.15 

GeV 2 



fit are shown in Tabic Q] and [Til together with the errors 
estimated by MINUIT. 

As the two different sets of Belle data are based on a 
different analysis of the same experimental events, they 
are strongly correlated. Therefore, we have treated them 
separately in our combined analysis of the HERMES, 
COMPASS and Belle data; the best fit values of Table 
H] are obtained by fitting the SIDIS results together with 
the Belle data on the cos(</?i +1P2) dependence, Eq. (|4"Tj) . 
while the values in Table |TT] originate from the Belle data 
on the cos(2</>i) dependence, Eq. (|59)) . We notice that the 
two sets of resulting best fit parameters are in full agree- 
ment within the uncertainties; this gives a good check 
of the consistency of the measurements and the stabil- 
ity of our analysis. In the sequel we shall present results 
and predictions based on the values of Table |T| the cor- 
responding results based on the values of Table [H] arc 
hardly distinguishable (examples of this arc shown ex- 
plicitely in Fig. [6] and in Fig. [8j right panel) . 

Our best fits of the experimental data from HERMES, 
COMPASS and Belle arc shown in Figs. H and © re- 
spectively. The central curves correspond to the central 
values of the parameters in Tabled while the shaded ar- 
eas correspond to one-sigma deviation at 90% CL and 
are calculated using the errors and the parameter cor- 
relation matrix generated by MINUIT, minimizing and 
maximizing the function under consideration, in a 9- 
dimensional parameter space hyper-volume correspond- 
ing to one-sigma deviation. 

The transversity distribution functions Ar?t(i,fci) 
and Axd(x, k±) as resulting from our best fit - Eqs. 1]13[) - 
(fTT|) and Tablc|T]- arc plotted as a function of x and k± in 



Fig. El for comparison, the Softer bound of Eq. (fl8|) is also 
shown, as a bold line. The solid central line corresponds 
to the central values in Table [I] and the shaded area cor- 
responds to the uncertainty in the parameter values, as 
explained above. 

Similarly, the resulting Collins functions 
A N Df av (z 7 p±) and A N D un f(z,p±) are plotted as a 
function of z - integrated over d 2 pj_, Eq. (|49|). and nor- 
malized to twice the unpolarizcd fragmentation functions 
- and as a function of p±_ in Fig. [HJ for comparison, we 
also show the Collins functions from Refs. [2J| and [2a ]. 
respectively as dashed and dotted lines (left panels), and 
the corresponding positivity bound (TR*)[) . The dashed 
lines in the right panels show the results corresponding 
to the parameters of Table [H] 
A few comments are in order. 

• In Fig. [7] we show the extracted transversity dis- 
tribution for u and d quarks. The a; dependence 
is based on the simple parameterization assumed 
in Eqs. (fil}j) and (TPo]) , which contain Nq ,a and 
(3 as free parameters; our result represents the 
first extraction ever of the transversity distribu- 
tions Atu(x) and Axd(x). 

• The k± dependence has been assumed to be the 
same as for the unpolarized distributions. The fla- 
vor dependence is contained in the coefficients Nq 
and in the proportionality of Axq(x) to [q(x) + 
Aq(x)]/2 = g+(a;), the number density of quarks 
with positive helicity inside a positive hclicity pro- 
ton. 

• Our results show that the transversity distribution 
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FIG. 4: HERMES experimental data @, [j| on the azimuthal 
asymmetry .4™HPs+pw for tt^ 1 production are compared to 
the curves obtained from Eq. (|20|) with the parameteriza- 
tions of Eqs. (|13|1 - (|17[1 . and the parameter values, determined 
through our global best fit, given in Table U The shaded area 
corresponds to the theoretical uncertainty on the parameters, 
as explained in the text. 



is positive for u quarks and negative for d quarks; 
the magnitude of Atu is larger than that of A-pd, 
while they are both significantly smaller than the 
corresponding Soffer bound. 
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FIG. 5: The measurements of A^ (0s+0h) , for the produc- 
tion of positively and negatively charged hadrons, from the 
COMPASS experiment operating on a deuterium target [l(| 
are compared to the curves obtained from Eq. (|20|l with the 
parameterizations of Eqs. (|13p ~ (|17p . and the parameter val- 
ues, determined through our global best fit, given in Table U 
The shaded area corresponds to the theoretical uncertainty 
on the parameters, as explained in the text. Notice the extra 
7r phase in addition to <j>s + <j>h in the figure label, to keep 
into account the different choice of the Collins angle, with re- 
spect to Trento [TJ and HERMES conventions, adopted by 
COMPASS Collaboration. 



The shaded regions in Fig. [7] show that both 
ATu(x,k±) and Ard(x, k±) are, considering the 
limited amount of data, already well determined. 
It is worth noticing that while the HERMES data 
alone tightly constrain the transversity distribution 
of u quarks, the addition of COMPASS data to the 
fit allows to better constrain the transversity dis- 
tribution function of d quarks. We have checked 
that fitting only HERMES and Belle data, ignor- 
ing the COMPASS results, still leads to a similar 
good x 2 /d.o.f.; the resulting functions would give a 
slightly worse description - when compared to the 
global fit - of the x dependence of A^s+M, as 
measured by COMPASS. This is mainly related to 
a less stringent determination of Axd(x, fcj_) in ab- 
sence of deuteron target data. Although their mea- 
sured azimuthal asymmetry is very small, the in- 
clusion of COMPASS data significantly contributes 
to the extraction of the transversity distributions. 
Different fitting procedures were earlier attempted, 
for example by fixing Axq = Aq or Axq = q\ [2q |: 
they lead to a slightly worse description of Belle 
data. 



• The extracted Collins functions are shown in Fig. 
[5] they agree with similar extractions previously 
obtained in the literature [H, [25[ . The shaded ar- 
eas indicate well constrained Collins functions for u 
and d quarks in the large (valence) z region, much 
smaller than their corresponding positivity bound. 



• We note once more that, in analyzing SIDIS data, 
we have neglected the contributions of the sea 
quarks and antiquarks (assuming the correspond- 
ing transversity distributions in a proton to vanish), 
taking into account only u and d flavors. In ana- 
lyzing Belle data and introducing the favored and 
unfavored Collins fragmentation functions, we have 
considered the contributions of u, d and s quarks, 
all abundantly produced in the e + e~ annihilation 
at Vs ~ 10 GeV. 



The partonic distribution and fragmentation func- 
tions are taken from Refs. [ij], [U| an d [3- The 
QCD evolution is taken into account in the unpo- 
larized distributions, in the unpolarized fragmen- 
tation functions and, following Ref. [13], for the 
transversity distributions. 
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FIG. 6: The experimental data on two different azimuthal cor- 
relations in unpolarized e + e~ — > h\hiX processes, as mea- 
sured by Belle Collaboration Q, are compared to the curves 
obtained from Eqs. (|44)l [A12] and (|60)l [Ao] with the param- 
eterizations of Eqs. (|14|l . (|16|l and (|17|l . The solid lines corre- 
spond to the parameters given in Table HI obtained by fitting 
the A12 asymmetry; the shaded area corresponds to the the- 
oretical uncertainty on these parameters, as explained in the 
text. The dashed lines correspond to the parameters given in 
Table [TT] obtained by fitting the Ao asymmetry. The agree- 
ment between the results obtained from the two fits shows 
the consistency between the two sets of Belle data and the 
solidity of our analysis. 
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FIG. 7: The transversity distribution functions for u and d 
quarks as determined through our global best fit. In the left 
panel, x Atu(x) (upper plot) and x AtcI(x) (lower plot), see 
Eq. ([5}, are shown as functions of x and Q 2 = 2.4 GeV 2 . The 
Soffer bound [23| is also shown for comparison (bold blue line). 
In the right panel we present the unintegrated transversity 
distributions, x Aru(x,k±) (upper plot) and x ATd(x,k±) 
(lower plot), as defined in Eq. (|13p . as functions of k± at 
a fixed value of x. Notice that this k± dependence is not ob- 
tained from the fit, but it has been chosen to be the same 
as that of the unpolarized distribution functions: we plot it 
in order to show its uncertainty (shaded area), due to the 
uncertainty in the determination of the free parameters. 
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FIG. 8: Favored and unfavored Collins fragmentation func- 
tions as determined through our global best fit. In the left 
panel we show the z dependence of the p± integrated Collins 
functions defined in Eq. (|38p and normalized to twice the cor- 
responding unpolarized fragmentation functions; we compare 
them to the results of Refs. [24| (dashed line) and [25| (dot- 
ted line). In the right panel we show the p± dependence of 
the Collins functions defined in Eq. (|14[) , at a fixed value of 
z. The Q 2 value is 2.4 GeV 2 , having assumed that the Q 2 
evolution of A N D is the same as that of D. The solid lines 
show the results based on the parameters of TableH] while the 
dashed ones show the results corresponding to the parameters 
of Table HTl In all cases we also show the positivity bound (|19[) 
(upper lines). 
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Finally, we explicitly list, for clarity and completeness, 
the kinematical cuts we have imposed in numerical inte- 
grations, according to the setup of the HERMES exper- 
iment: 

0.2 < ^ < 0.7, 0.023 < x B < 0.4 , 

0.1 < y < 0.85, Q 2 >lGeV 2 , (61) 

W 2 > 10 GeV 2 , 2 < E h < 15 GeV , 

the COMPASS experiment: 

0.2<2h<l, 0.1 < y < 0.9, 

Q 2 > 1 GeV 2 , W 2 > 25 GeV 2 , (62) 

and the Belle experiment 

- 0.6 < cos 0i ab < 0.9 , Q T < 3.5 GeV , (63) 

where Ouh is the polar production angle in the laboratory 
frame (related to the scattering angles 8 and #2 used in 
this paper) and Qt is the transverse momentum of the 
virtual photon from the e + e _ annihilation in the rest 
frame of the hadron pair [23| . 
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FIG. 9: Predictions for the single spin asymmetry A s ^ s+M 
as it will be measured by the COMPASS experiment operating 
with a transversely polarized hydrogen target. For the extra 
7r phase in the figure label see the caption of Fig. 5. 



For a neutron target at 6 GeV the cuts read: 



V. PREDICTIONS FOR ONGOING AND 
FUTURE EXPERIMENTS 

We can now use the transvcrsity distributions and the 
Collins functions we have obtained from fitting the avail- 
able HERMES, COMPASS and Belle data, see TableU to 
give predictions for new measurements planned by COM- 
PASS and JLab Collaborations. 

The transverse single spin asymmetry A s ^^ s+( ^ h ^ will 
be measured by the COMPASS experiment operating 
with a polarized hydrogen target (rather than a deu- 
terium one). In Fig. [9] we show our predictions, obtained 
by adopting the same experimental cuts which were used 
for the deuterium target, see Eq. ([62]) . Notice that this 
asymmetry is found to be sizeable, up to 5% in size. 

The JLab experiments will measure j\ s ^i<Ps+4>h) f or 
pion production off transversely polarized proton and 
neutron targets, at incident beam energies of either 6 
or 12 GeV. The kinematical region spanned by these ex- 
periments is very interesting, as it will enable to explore 
the behavior of the transversity distribution function at 
large values of x, up to x ~ 0.6. The adopted experimen- 
tal cuts for JLab operating on a proton target at 6 GeV 
are the following 



0.4 < z h < 0.7, 
0.1 < x B < 0.6, 
Q 2 > 1 GeV 2 , 
1 < E h < 4 GeV , 

whereas for a beam energy of 12 GeV they are 



0.02 < P T < 1 GeV, 
0.4 < y < 0.85, 
W 2 > 4 GeV 2 , 



(64) 



0.4 < z h < 0.7 , 
0.05 < x B < 0.7, 
Q 2 > 1 GeV 2 , 
1 < E h < 7 GeV . 



0.02 < P T < 1.4 GeV, 
0.2 < y < 0.85, 
W 2 > 4 GeV 2 , 



(65) 



0.46 < z h < 0.59, 

0.68 < y < 0.86, 

5.4 < W 2 < 9.3 GeV 2 , 



0.13 < x B < 0.40, 
1-3 < Q 2 < 3.1 GeV 2 , 
2.385 < E h < 2.404 GeV 



(66) 

whereas for an incident beam energy of 12 GeV they are: 



0.3 < z h < 0.7, 
0.34 < y < 0.9, 
W 2 > 2.3 GeV 2 . 



0.05 <x B < 0.55, 

Q 2 > 1 GeV 2 , (67) 



Our corresponding predictions, according to Eq. (|20|) 
and our extracted transvcrsity and Collins functions, arc 
shown in Figs. [TU] and [TT1 

It is important to stress that, as the large x re- 
gion is not covered by the HERMES and COMPASS 
experiments, our predictions for the x dependence of 
^sin(0 s +0h) arc vcr y sensitive to the few available data 
points from HERMES and COMPASS at moderately 
large x values. As a consequence, the predictions for 
the JLab experiments may vary drastically in the region 
0.4 < x B < 0.6, as indicated by the large shaded area 
in Figs. [TO] and [TTJ On the contrary, the results on the 
Pt and z^ dependences are more stable, as they only de- 
pend on the transversity distribution function integrated 
over x. 

Finally, we compute the azimuthal asymmetry 
j^sin(<ps+4>h) f or the production of K mesons and com- 
pare it with existing HERMES results @,[9[. These data 
have not been included in our best fit, as they might in- 
volve the transversity distribution of strange quarks in 
the nucleon, which we have neglected for SIDIS data on 
7r production. We show our results in Fig. [TOJ obtained 
using the extracted u and d transversity distributions. 
Again, we have used favored (A^fl^+^t) and unfa- 
vored (A N D K - / u j , A N D K ± / d r ) Collins functions, as in 
Eqs. (fT4"|) . (fTO[) and (fTT)) . For these we have used the same 
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10: Predictions for the single spin asymmetry 



^sin(0 s +</) h ) ag y. w jjj k e measured at JLab operating on polar- 
ized hydrogen (proton, upper plot) and He 3 (neutron, lower 
plot) targets at a beam energy of 6 GeV. 
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FIG. 11: Predictions for the single spin asymmetry 
^sin(0 s +^ h ) ag y. w jjj k e measured at JLab operating on polar- 
ized hydrogen (proton, upper plot) and He 3 (neutron, lower 
plot) targets at a beam energy of 12 GeV. 



parameters Nq 



7, 8 and M of Table [Q with the appro- 
priate unpolarized fragmentation functions D K ± j q [19( . 

We notice that our computations are in fair agreement 
with data concerning the K + production, which is pre- 
sumably dominated by u quarks; instead, there seem to 
be discrepancies for the K~ asymmetry, for which the 
role of s quarks might be relevant. New data on the 
azimuthal asymmetry for K production, possible from 
COMPASS and JLab experiments, might be very helpful 
in sorting out the eventual importance of the sea quark 
transversity distributions in a nucleon. 



VI. 



COMMENTS AND CONCLUSIONS 



We have performed a combined analysis of all exper- 
imental data on spin azimuthal asymmetries which in- 
volve the transversity distributions of u and d quarks 
and the Collins fragmentation functions, classified as fa- 
vored (when the fragmenting quark is a valence quark 
for the final hadron) and unfavored (when the fragment- 
ing quark is not a valence quark for the final hadron). 
We have fixed the total number of 9 parameters by best 
fitting the HERMES, COMPASS and Belle data. 
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FIG. 12: Our results, based on the extracted transversity and 
Collins functions, for the azimuthal asymmetry A B ^^ s+ ^ h ' 
for K ± production, compared with the HERMES experimen- 
tal data HQ. 



All data can be accurately described, leading to the 
extraction of the favored and unfavored Collins functions, 
in agreement with similar results previously obtained in 
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the literature [25[. In addition, we have obtained, 
for the first time, an extraction of the so far unknown 
transvcrsity distributions for u and d quarks, Atu(x) 
and AtoI(x). They turn out to be opposite in sign, with 
|Ard(x)| smaller than |Aru(x)|, and both smaller than 
their Soffer bound [IJ. 

The knowledge of the transversity distributions and 
the Collins fragmentation functions allows to compute 
the azimuthal asymmetry A^Ote+^s) for any SIDIS pro- 
cess; we have then presented several predictions for in- 
coming measurements from COMPASS and JLab exper- 
iments. They will provide further important tests of our 
complete understanding of the partonic properties which 
are at the origin of SSA. Data on K production will help 
in disentangling the role of sea quarks. 

Further expected data from Belle will allow to study in 
detail not only the z dependence of the Collins functions, 
but also their p± dependence. The combination of data 
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